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Q_|. Abstract 

1) ; 

. I Two-photon decays probe the structure of mesons and represent an important con- 

tribution to hadronic light-by-light scattering. For the neutral pion, the decay amplitude 
tests the effects of the chiral anomaly; for a heavy quarkonium state, it measures the 
magnitude of its wavefunction at the origin. We rederive the expression of the decay 
amplitude in terms of a Euclidean correlation function starting from the theory de- 
fined on the torus. The derivation shows that for timelike photons the approach to the 
' infinite-volume decay amplitude is exponential in the periodic box size. 
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$_i , 1 Introduction 



Photon-hadron interactions constitute a venerable subject that continues to be impor- 
tant in 21st century particle physics. Historically, the inner structure of the proton was 
investigated in elastic scattering experiments of electromagnetic probes pQ. Roughly 
a decade later, the study of deeply inelastic processes [2j[3] led to the development of 
the parton model, thereby making a decisive contribution to the advent of QCD as the 
theory of the strong interactions. At the same time, the decay ttq — > 77 led to the 
realization that the flavor-singlet chiral symmetry of the classical field theory is broken 
by quantum effects [HE]- in this article we rederive the expression for the hadronic 
matrix element that determines the two-photon width of the neutral pion, starting from 
QCD defined on a finite torus. This is the context in which almost all lattice QCD cal- 
culations are performed. An important conclusion of this analysis is that the finite-size 
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effects on the type of hadronic matrix elements that several lattice QCD collaborations 
are computing [5HH] are exponentially suppressed. 

We assume that QCD has exact isospin symmetry prior to coupling hadrons to 
photons. We consider a neutral, pseudoscalar meson of mass M which is stable in 
QCD and which we refer to as the 'pion', although we need not assume that it is a 
pseudo-Goldstone boson. The basic idea is to think of the square module of the decay 
amplitude as the pion pole contribution to the light-by-light scattering amplitude. We 
will be led to consider a center-of-mass energy off by order e 2 from the pion pole (e is the 
electromagnetic coupling constant). At this energy, the 77 cross-section is of order e 2 
rather than e 4 . For the sake of the argument, we will consider massive vector bosons v 
with a finite mass M v < M/2 and denote the corresponding field theory by mQED. The 
QCD part of the ttq decay amplitude into two timelike photons is reinterpreted as the 
QCD part of the ttq decay amplitude into two massive, on-shell vector bosons. We can 
then use the general formalism developed by Liischer [9fl0] to establish a correspondence 
between the (discrete) spectrum of mQED+QCD on the torus and the vv scattering 
phase. The ttq — vv interaction leads to a splitting of the ttq mass and a two-vector- 
boson energy level, given by matrix elements of the finite-volume Hamiltonian. The 
relevant matrix element can then be rewritten in terms of a time-ordered product of 
vector currents in the finite- volume theory. Since the relation between the finite- volume 
spectrum and the scattering amplitude holds up to exponentially suppressed corrections, 
we find that this conclusion also holds for the module of the ttq — > 7*7* amplitude. 

Once the finite-volume corrections have been treated in this way for timelike pho- 
tons, it is a seemingly innocuous step to extend the conclusion about the finite-size 
effects to the case of real (lightlike) photons. Rather than giving a formal proof, we 
check that the relation between the finite- volume matrix element and the infinite- volume 
matrix element holds even for massless photons in two simples field theories. 

One might at first think that, photons being massless, the finite- volume effects will 
only be suppressed by a power of the linear extent of the torus. However, at leading non- 
trivial order in the electromagnetic coupling, the photons are merely external particles in 
the reaction of interest. Their interaction is of purely hadronic origin, and has a range 
of order the inverse pion Compton wavelength. Inelastic channels are kinematically 
allowed; but they are suppressed by the smallness of the electromagnetic coupling. We 
note that the partial wave expansion which underlies Liischer's formalism is also valid 
for relativistic and even massless particles (see for instance section 3.7). 

Our treatment of the problem bears a strong similarity with the formulation of the 
K — > 7T7T decay on the torus by Lellouch and Liischer [12] . In some ways, the ttq — > 77 
amplitude is simpler: there are no final-state interactions at order e 2 . In other ways, 
it is more complicated: the decay amplitude is of second order in the photon-hadron 
interaction Hamiltonian. The latter induces a shift in the pion mass which is also 
formally of order e 2 , and which is different in finite and infinite- volume for massless 
photons (power-law corrections are expected). A second minor complication arises due 
to the two possible helicities of the final-state photons (either both positive or both 
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negative), since there are in principle two open elastic channels. However, due to parity 
conservation, the two channels are related by a phase difference and one can reduce the 
problem to a single-channel problem. 

The structure of this paper is as follows. The master equation (I2,20|) relating 
the finite-volume matrix element to an infinite-volume matrix element is derived for 
timelike photons in section 2, where it is also brought into a form suited for lattice QCD 
calculations. In section 3, the pion decay amplitude is calculated in two simple field 
theories coupled to massless photons, one being an effective description for light quarks 
and the other for heavy quarks. In section 4, the corresponding finite- volume matrix 
elements are calculated directly, thus providing a check of the master equation. Unless 
otherwise stated, we use the notation and conventions of Peskin and Schroeder [13], in 
particular the metric convention g^ v = diag(l, —1,-1,-1). 

2 General considerations 

We introduce the amplitude for the decay 7To — > 7*7*. In standard notation (see for 
instance jT3] section 5.1), the amplitude is given by e 2 times the QCD matrix element 

A^(p,q) = i j 'd 4 xe^(0|T{i /1 (x)^(0)}Ko(p)) re i (2.1) 
= e^ a(3 q a pP T(ml,q 2 ,(p-q) 2 ), (2.2) 

where is the electromagnetic current of hadrons. We have made explicit the use of 
the 'relativistic' normalization of states, 

rc\(MP )\MP)) rei = 2E P (2tt) 3 5^ (p - ft) . (2.3) 

The quantity J- is the off-shell ttq —> 77 form factor. We will only consider the situation 
where the pion is on-shell up to higher order corrections and (for simplicity) the case of 
two photons of same invariant mass, 

A = ^*(qy a *(p-q) I d 4 xe^(0\T{j,(x)j u (0)}\Mp))re h a = +, (2.4) 

with p 2 = M 2 and q 2 = (p—q) 2 - We generically denote by e a ^(k) the polarization vector 
of the photon. Only the a = ±1 polarizations will play a role, because the longitudinal 
polarization a = does not contribute to the amplitude. The pole contribution to the 
Lorentz-invariant light-by-light scattering amplitude is then 

iM( 7 +7+ -> ^0 -> 7+7+) = ~e A A* p2 % _ M2 A. (2.5) 

Except in the propagator, the four-momentum of the pion can be considered to be 
on-shell. Due to angular momentum conservation, the two photons must have the 
same helicity. Due to parity conservation, there is a relative minus sign between the 
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amplitudes for two positive-helicity photons and the amplitude for two negative-helicity 
photons in the final state. This information is already encoded in the parametrization 
(12~2D . 

We will need the lowest term of the partial-wave expansion of the light-by-light 
scattering matrix. Due to the two possible helicity states of the two photons in the final 
state, we are really dealing with a (special case of the) two-channel problem. In the 
s-wave, the partial wave expansion of the amplitude (|2.5h reads (see appendix A) 

M = \- lGirE ■ t , to = ^ (e 2iSo - l) • (2.6) 

2.1 Discrete energy eigenstates on the torus 

Our goal is to extract the decay rate of the neutral pion into two photons, ir° — > 77, 
from stationary physics on the torus. For the sake of the argument we will however 
consider massive vector bosons v with a finite mass M v < M/2. The QCD part of the 
7To decay amplitude into two timelike photons is reinterpreted as the QCD part of the 
7To decay amplitude into two massive, on-shell vector bosons. 

The linear size L of the torus will be assumed to be in the range 



2vr 

T 



A ms- ( 2 - 7 ) 



In current state-of-the-art lattice QCD calculations, L is typically 3 to 5 fm, in which 
case 2tt/L = 1240MeV/L[fm] is in the ballpark of (|2.7p . Also, the pion is assumed to 
fit well into the box, so that the mass difference between the finite-volume theory and 
the infinite- volume theory can be neglected |15j . 

We will restrict our attention to the center-of-mass system throughout this paper. 
We work perturbatively in the electromagnetic coupling e. Initially we set e = 0, so 
that the vector-boson states and the hadronic states are decoupled. The spectrum of 
single-vector-boson states is 



jw fc k = ^n, n G Z 3 | , oo k = ^k 2 + M$, 



(2i 



and the spectrum of two-photon states with vanishing total momentum is 

2u k k = jn, n G Z 3 | . (2.9) 
The interaction Hamiltonian reads 

V = +e j dxj»(x)Ap(x). (2.10) 

The energy levels will be affected by the electromagnetic coupling. Consider first non- 
degenerate states. At first order in ordinary time-independent perturbation theory, the 
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non-degenerate energy levels will change by AEn^ = (n\V\n) . We normalize the finite- 
volume states to have unit norm, (n\n) = 1. However, the interaction Hamiltonian 
changes the vector-boson number by exactly one unit, therefore the matrix elements 
between states containing the same, definite number of vector bosons vanish. At second 
order, 

If the unperturbed state \n) is a non-degenerate, purely hadronic state, then the states 
\m) that contribute non-trivially to the sum in Eq. (|2.1ip are states that contain one 
vector boson and have an arbitrary hadronic content. For a state \n) containing N n 
vector bosons, contributions come from states with N n ± 1 vector bosons. 

Here we will be interested in the case where the box size is tuned so that the 
neutral pion mass coincides with the energy of exactly one two-vector-boson state with 
vanishing total momentum. As is well known, the corrections to the spectrum are found 
by diagonalizing the perturbing Hamiltonian in the degenerate subspace, which amounts 
to diagonalizing a finite matrix V nn i . In this case, the 2x2 matrix V nn i vanishes since 
there is no pair of states for which the two states differ by a vector boson number of 
one. To find the leading effect of e.m. interactions on the degenerate |7ro), \vv) states 
we therefore have to go to second order of degenerate perturbation theory. The recipe 
is then to diagonalize the matrix ([16]. paragraph 39) 

where the sum over m extends over all states not in the degenerate subspace. Since 
V nn > vanishes, the matrix W nn > is 0(e 2 ). We label the neutral pion at rest by '1' and 
the two-vector-boson state by '2'. For instance, the matrix element W22 corresponds 
to the transition of one vector boson to a hadronic state and back, while the other one 
remains as a spectator, a vacuum polarization effect due to the hadronfl The matrix 
element Wu correspond to a mass correction for the pion due to the weakly coupled 
vector bosons. Its dependence on the volume is exponentially suppressed^]. The matrix 
W nn i takes the form 

and the final energy eigenvalues are thus predicted to be, to 0(e 2 ), 



E± = M + AE± = M + U L M + 5 L E V ± \ (U L M - 5 L E V ) 2 + \W 12 \ 2 . (2.14) 



lr There is a further potential effect, stemming from the fact that the electromagnetic interaction 
breaks isospin symmetry, which is the mixing of the no with the rj meson and other, not necessarily 
isovector states. Since the no — r\ mixing starts at 0(e 4 ), we can neglect it. 

2 For instance, for a non-relativistic bound state described by a wavefunction ^(r), it would be 
AM = —e 2 J d 3 r \tp(r)\ 2 e in infinite volume and would only differ by an exponentially small 

amount on the torus, provided the bound state fits well into the box. 
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2.2 Light-by-light scattering near the pole 

Hadronic light-by- light scattering is generically of order e 4 , however for energies a dis- 
tance e 2 from the pion pole, the scattering amplitude and the phase shift are of order 
e 2 . The decay width of the tto is of order e 4 and therefore negligible in the following 
considerations. However, we must take into account the 0(e 2 ) shifts in the ttq mass and 
the photon mass. 

We treat 5q as a function of the vector boson momentum. For a small phase shift, 
Eq. ([23]) and QDft yield 

fce 4 |.4| 2 

6o{k±) = - 16^(AE ± - 6 00 MY (2 - 15) 

2.3 Relation of the discrete energy levels to the ttq — > vv amplitude 

In the finite- volume theory, we work in the basis formed by the state 

im-> = ^=E £ ™L°LJ°> ( 2 - 16 ) 

v n feen n cr=±i 

and the state for which the relative sign between the two helicity terms is 

positive. The set of momenta appearing in Eq. (|2.16p is 

n n = {k = 2irz/L | z G Z 3 , z 2 = n}, 

where we restrict ourselves to a value of n < 6 in order to avoid multiple degeneracies, 
and v n = |f2 n | is the cardinality of that set. Only couples to the ttq, so that 

we can ignore the state \ (vv)+) to the order e 2 we are working at. The discrete set of 
momenta T comprises all the images under the cubic group of one particular allowed 
momentum on the torus. 

Using Liischer's finite- volume formalism [9,10j, a change in the effective momentum 



k defined via the two- vector-boson spectrum, E = 2\J M% + k 2 , from the free-field values 
(|2.9p corresponds to an s-wave scattering phase [T2"1[TT] 

A«t)=-(rf(,) + t -|i) T , , = - (2.17) 

with (j) a tabulated kinematic function. Here 5o(k) vanishes before the vector bosons 
are coupled to hadrons. Note that the electromagnetic shift in the ttq pole is contained 
in A<5o- We use the derivative of (ft at a point where q = \z\ for some vector z G Z 3 , 
q<j)'{q) = (2vr) 2 ^. One thus finds 

5o(k±) = -^Ak ± . (2.18) 
2irv n 
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Now the change in the effective momentum of the two vector bosons is related to the 
spectrum via 

M 

Ak± = — (AE± - 25 L E V ) . (2.19) 
4/c 

Combining Eq. (|2,15p . (|2. 18|) and (|2.19p . we obtain the master relation 



9 T 3 /\#3 or3 71^3 

\A\ 2 = — j— ■ (AE ± - 6 OQ M)(AE ± - 25 L E V ) = — j— • \W 12 \ 2 - (2.20) 



We have used the fact that 5lM = 5ooM up to exponential corrections, and neglected 
the latter. 



2.4 Explicit expression for the matrix element Wyi 

Using the spectral representation, one verifies that Wi%, defined in Eq. (|2.12p . is given 

by 

Wn = -\J" J* <M-|T{U(W0)}|£>, (2.21) 

as long as there is no contributing intermediate hadronic state with energy E^ < M/2. 
The Euclidean time evolution, V(t) = e m V(0)e~ Ht , is dictated by the full QCD+mQED 
Hamiltonian H. The Wick contractions of the vector bosons can be carried out straight- 
forwardly by noting that in equation (|2.12p . the intermediate states contributing are di- 
rect products of an arbitrary hadronic state with a one-vector-boson state and by using 
equation (|2.16p and the plane-wave expansion of the gauge field. One then reduces W±2 
to a pure QCD matrix element, 

/ p u) k \xo\+ik-x 
d A x — (0\T{f(x)f(0)}\B). (2.22) 

Taking into account the relativistic normalization of states (s/2ML' i \B) — > |i?) re i.) and 
inserting (|2.22[) into (|2.20p shows that, possibly up to an overall phase (p, the ttq — > vv 
amplitude is given by the analytic continuation of Eq. (|2,4p . 

A = e l *e* +ll (-k)e%(k) J d 4 x e "*l*°l+** (0\T{f(x)f(0)}\B) iel . (2.23) 

This expression is suitable for an implementation in lattice QCD [18], see [6HH] for recent 
calculations. 



3 The pion decay amplitude in two simple field theories 

In the previous section, relation (I2.20p was derived for massive vector bosons and then 
reduced to a pure QCD matrix element, Eq. (|2.23p . In the latter equation, the only 
memory of the dispersion relation of the vector bosons is the relation between Uk and 
the spatial momentum k. Since the vector bosons merely appear as external legs in 
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the amplitude, we expect that the volume corrections affecting the calculation of the 
amplitude A in finite volume based on Eq. (|2.23p are exponentially small in the box 
size, even for the case of real photons where u>k = 

In order to check this expectation, in this section we calculate, in two tractable 
field theories, the amplitude A for on-shell photons in the final state. We then use 
relation (|2.20p to predict the finite- volume quantity |VFi2| 2 - In the next section, the 
same quantity will be computed directly in finite volume as a check of Eq. (|2.20p . 



3.1 Theory I: chiral Lagrangian 

In chiral effective theory the decay of a neutral pion of mass M into two photons is 
driven by the effective Lagrangian (see for instance |19j . chapter 22) 

C a = 9 vr e^F^F pa . (3.1) 

In the chiral limit the Abelian chiral anomaly [H[5] predicts the value of the coupling [20] , 

N c e 2 

9 = 48^; (3 - 2) 

with F n ~ 184MeV the pion decay constant and N c the number of colors, however we 
will not assume any particular value for g. In this theory the module of the invariant 
amplitudes for ttq decay then amounts to 

\M(B -> 7+7+ )| 2 = \M(B -> 7_ 7 _)| 2 = e A \A\ 2 = 16<? 2 M 4 , (3.3) 

which leads to a full 77 width of T = M s g 2 /ir. 

In order to verify (|2.20p . we thus have to show by a direct calculation that the 
module of the finite-volume matrix element W12 is given by 

1^1' (3.4) 

3.2 Theory II: a non-relativistic bound state of two massive fermions 

Consider first the amplitude for the fermion-pair annihilation // — > 77. The fermions 
have a mass m. The initial momenta of the fermions are p and p' and the final momenta 
of the photons are k and k' . There are two diagrams at treelevel, yielding 



,2 



M = -e^kfe^k')* v(p') 



2 



7 — + 



u{p) (3.5) 



p ■ k p ■ k' 

We have used p 2 = m 2 and k 2 = and the spinor identities 

(J> + m)j u u(p) = 2p u u(p), u{p)^f v (^> + m) = u(p)2p u '. (3.6) 
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Now we choose the center-of-mass frame and consider the limit of non-relativistic 
fermions. Using standard spinor technology one obtains (for an initial state with one 
spin-up fermion and one spin-down fermion) 

He 2 

M = — (kxe a (k)*)-e a ,(-k)*. (3.7) 

m 

For instance, if k = \k\e3 and e a (k) = ^=(l,ia, 0), e CT /(— fc) = -^(1, — zcr', 0), then 

i(c 3 x e a (k)*) ■ e a/ (k')* = -aS aal . (3.8) 
Following ([13], chapter 5), we describe a pseudoscalar bound state a^| 

The amplitude for B — > 77 is then 

A4(B-» 7( /v) = V2M I ^2-^M.{M{p^ _p^ ry ) + ^(p4. > -pt-^77)) 

y (Z7TJ J 2m 

= -a<W' ^£ *(* = 0). (3.11) 



As a check, the decay rate of the bound state of mass M into two photons now reads 
' d " k f ''' k ' - {2^ A 5^\ P + p'-k-k')\M{B^ 11 )\ 2 



2M J (2vr) 3 2|fc| J (2vr) 3 2|fc 



= ^|* (a; = 0)| 2 . (3.12) 

In r, we have implicitly included a factor 2 stemming from the sum over final polar- 
ization states of the photons, J2 a J2 a ' ^ = 2> but also a factor h to remove indistin- 
guishable final-state configurations of the photons. For muonium in the ground state, 
we have ^S>(x = 0) = (/x 3//2 e 3 )/(87r 2 ), with fi = m/2 the reduced mass, and the width 
becomes the well-known result T = a 5 m/2. 

By comparison with Eq. (|2.4p . the result (|3.1ip implies in particular 

\A\ 2 = -\y(x = 0)\ 2 . (3.13) 
m 

In order to verify (|2.20p . we thus have to show by a direct calculation that the module 
of the finite-volume matrix element W\2 is given by 

W 2 = ^(-«)| 2 (3-14) 



3 The first factor in H3.10|) converts the state from a non-relativistic normalization to a relativistic 
one and we have defined 

\pt -pi) = <(V-p(-] 1 0), I p i, -p t ) = <(V^ ( +) I °)> ( 3 - 9 ) 

where a F (a F ) is the annihilation operator of the particle (antiparticle). The relative sign between the 
two terms in (|3.10[) can be checked by calculating the (non- vanishing) amplitude (O|'075'0|S). 



9 



4 Analytic calculation of finite-volume matrix elements 



In order to check the master formula (|2.20p , we need to determine by a direct calculation 
the norm of the matrix element W\2- We work in the Euclidean theory and adopt the 
corresponding notation. In particular, the Euclidean Dirac matrices are related to the 
Minkowski space matrices by 7q = 7 and 7^ = — i~/ k , 7^ =75, but we will not display 
the superscript 'E' explicitly. 



General considerations and the spectral representation 

In order to compute W12, we consider the gauge-invariant correlation function 

G(t) = (V(t)O 5 (0)} . (4.1) 
A gauge-invariant choice of operators would be 

V ^ = J d 3 xe tlvpa F^(t,x)F pa (t,x) (4.2) 

= JJ/2- J d 3 x4,(t,x) 75 iP(t,x). (4.3) 

In order to simplify the perturbative calculation we will however make a slightly different 
choice in the next section. We will compute G{t) perturbatively in the electromagnetic 
interaction at order e 2 . The first non-trivial contribution to G(t) occurs at order e 2 , 

G(t) = Go(t)+e 2 G 2 (t) + e 4 G 4 (t) + ..., (4.4) 

(hit) = ±(nt)(fd 4 zV(z)) 2 O 5 (0)) Q , (4.5) 

V(x) = -ieA^x^ix^ipix), (4.6) 

V(t) = [ d 3 x V(t,x). (4.7) 



In these equations the (Euclidean) time dependence of the operators is dictated by the 
non-QED part of the Hamiltonian. 

On one hand, the correlation function can be interpreted in terms of the spectral 
representation. In order to isolate the desired matrix element, we notice that the con- 
tribution of interest comes about when the interaction Hamiltonian f d 3 x V is inserted 
between time and t, 

G 2 (t) = V (o\r(t) n) f ds e - En{t ~ s) / S ds' (4.8) 
n,n> X ' ' Jo Jo 

^2{n\V\l) e- E ^ s ~ s 'Hl\V\n') e"^' s ' In' O 5 {0) 0) + . . . 



We assume that the torus size has been tuned so that the states |n) = I77) and \n') = \B) 
are degenerate with an energy E = 2\k\. Carrying out the integrals over s and s', their 
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contribution to G 2 {t) will have the following dependence on t, 

( 71 \V\l) (l\V\B) 



G 2 (t) = te- Et (V\ 17 )(B\0 5 )[J2 



Ei — E 



+ ■■■ 



where for brevity we have defined 



(7>| 77 ) = (0 7>(0) 



O 5 (0) 



. 



(4.9) 



(4.10) 



{B\0 5 ) = (B 

By comparing to Eq. (j2. 12|) . we see that we can extract the matrix element W\ 2 from 
G 2 (t) = -W 12 ■ (V\jj) ■ (B\0 5 )te- Et + ... (4.11) 

The overlap factors 

We use creation and annihilation operators that are unit-normalized, [%>sj a L s /]± = 
5 ss '5 P pi where the plus sign is for bosons and the minus sign for fermions and s represents 
all discrete indices characterizing the particles. The unit-normalized energy eigenstates 
we are interested in are the following. The two-photon state is defined in Eq. (|2.16p . 
For Theory I and II respectively, the 'pseudoscalar meson' states are 



\B) 
\B) 



4=o 1°)' 
1 



E A , s/ Ft Ft .Ft 



Ff Ff 
-)°-p(+) 



(4.12) 
(4.13) 



We will use the following spatial Fourier representation and normalization for the wave 
function, 



• -l^V/'-'ui,,/,). 



d s x \iff(x) 



L 3 



1. 



p p 
Our interpolating operator for the two-photon state will be 

V = Ai(—k)Aj(k), with i = 1, j = 2, k = |fe|e 3 . 

Using the plane- wave expansion of the gauge fields, 



one finds 



(^|(77)-) 



iL 3 



2u k Jv,, 



(4.14) 



(4.15) 



(4.16) 



(4.17) 
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for the matrix element. The matter fields have the expansion 

^o(x) = ^^==(a^ e ^ + a-t e -^) ; (4.18) 

^ = Y.^^Y.Usn s {p)e^ x + ap s v s (p)e~^), (4.19) 
P V ltj v L s 

^) = Y.^^Y.(<sV s (p)e^ x + a F p }u s {p)e-^) (4.20) 
p V V s 

with the normalization u r (p)u s (p) = 2m5 rs and v r (p)v s (p) = —2m5 rs . In theory I, the 
interpolating operator is simply the pion field, O5 = 7ro, and the overlap is 

; (4,21) 

in theory II, the interpolating operator is chosen to be O5 = — Y>(p)75i/>(— p) and the 
overlap is 

(5|0 5 ) = **(-p). (4.22) 

Perturbation theory in the time-momentum representation 

On the other hand, can be computed using the standard techniques of Feynman 
diagrams, using the time-momentum representation of the propagators. The fields are 
spatially Fourier transformed according to 

A^(t,k) = J d 3 x A^(t,x) e- ikx , 7r (t,fc) = J d 3 x ir (t,x) e~ ik ' x , (4.23) 

ip(t, fe) = y d 3 x ip(t, x) e~ ikx , tP(t,k) = j d 3 x tp(t,x) e~ ikx . (4.24) 

In this representation, the photon propagator in Feynman gauge reads 

-w k \t-t'\ 

(A^t, k) A u (t', -k')) = 5^L 3 5 kk , — , (4.25) 

and the matter fields have the propagators 



e -E p \t-t'\ 

e -E P \t-t'\ 



(n (t,p)*otf, -*>')) = L ^pp> ~ 2E , (4.26) 



{*/>(t,p) j(t', -p')) = L 3 5 PP , (E p sgn(t - t') 70 - ip ■ 7 + m) . (4.27) 

The fermion-photon vertex is —iej^ (fi corresponds to the Lorentz index of the photon 
and e is the charge of the fermion). 
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4.1 Theory I: finite- volume calculation 

At first, we will not make the assumption that the two-photon state is necessarily 
degenerate with the pseudoscalar matter state. The correlation function at order e 2 
with all fields written in the time-momentum representation reads 



G 2 (t) = ig^ upa {A x {t,-k)A 2 (t,k) (4.28) 
ds-i E^'^TI E F flu (s,l)F pa (s, -(I + <?)) vr (0,p = 0)' 



Q I 



After all Wick contractions have been performed, one arrives at 

G 2 (t) = _?|^L f°° ds sgn(t - s) e-^l.l+a-fclt-l). ( 4 .29) 
Now carrying out the integral over s yields 



ds sgn(t - s) e -(ms\+2u k \t-s\) = £ e e e_ 

B v ' 2u k + M 2oj k — M 



where the first term comes from the region s < and s > t and the second from 
< s < t. If we now take the limit ui k — > M/2, the first term vanishes and the second 
becomes te~ Mt . Thus the correlation function becomes, for k% = M/2, 

G 2 (t) = -^te- M K (4.31) 

Comparing it with Eq. (|4.11j) and dividing out the overlap factors (|4.17p and (|4.2ip , we 
indeed arrive at Eq. (I3.4p . thus confirming the master relation (|2.20p . 

4.2 Theory II: finite-volume calculation 

The correlation function at order e 2 with all fields written in the time-momentum rep- 
resentation reads 

e 2 / 

G2(t) = -^{A(t, -k) Aj(t, k) (4.32) 



/> 2? E f E Ms, e) jy± E h E ^ Q/) 



(75)^ 



The Wick contractions of the gauge fields lead to 

2 poo poo p -Wfc|t-s| p-w fe |t-a'| 1 1 

^(s, £)7V(«, + fc)) qh j ^(s', k-q) ^ a (0, -p)#(0, p)' 
+(i Hj, fe — > -k). 
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The fermion propagator joining the two insertions of the vector current is treated as a 
free propagator, 



oo e -u k \t-s\ e -u k \t-s'\ 



1 g-Sq+fc|s-s'l 



Z3 E * 2 ^| fc ' Sgn(s " S ' )7 ° + ^f(9 + k) + m) sy 

^ 7 (s, q)^' (s', -g) Va(0, -p)^(0,p) 
fc^ -fc). (4.35) 

We can now use expression (jB.6[) for the fermion four-point function that takes into 
account the presence of a bound state in the spectrum (see appendix B) and obtain 

G 2 (t) = 2e 2 ^(-p)* (u (+) (-p) 75 « ( " ) (p) + u (_) (-P)75« (+) (p)) (4.36) 

/s e -( J j k \t-s\ e -u> k \t-s'\ e -E B s' e -E q+k (s-s') e -E q (s-s') _ 
ds' — ^(-q) 
oo 2wfc 2wfc 2 • 2E p z — ' 2E q+ k 2E q 



+ fe -> -fc) + • • • 

The dots stand for terms which do not contribute to the contribution proportional to 
t e~ Et that we are interested in. Now since the wavefunction ^(q) is concentrated at 
momenta much less than m in norm, we can approximate all the factors at leading order 
by evaluating them at q = 0. After this operation, the first two lines in Eq. (14. 36ft are 
an even function of fc, hence we can evaluate the spinor matrix element where we now 
set i = 1, j = 2 and k = me^, 



J^^Wf^kSgnis - s')-f + i-y ■ k + mjj j u~ r {0) + {i ^ j, fc -fc) (4.37) 

r=± 

= -2im v r (0) t 7 5 u" r (0) = Urn 2 , 

r=± 

u (+) (0)7 5 v { " ) (0) + (0)7 5 ?j {+) (0) = -4m. (4.38) 
Now the integrals over s and s' can be carried out, yielding (for i = 1 and j = 2) 

e 2 L 3 te~ 2mt 

G 2 (t) = -i¥(0) vf(-p)* ^ + . . . (4.39) 

We assume that the wavefunction of the bound state on the torus only differs by a 
negligible amount from the infinite volume wavefunction. For L sufficiently large this 
assumption is justified by the general analysis [15]. Comparing with Eq. (|4.1ip and 
dividing by the overlap factors (|4.22p and (|4.17p . we arrive at the expected formula 
(I3.14p and have thus confirmed the master relation (I2.20p . 
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5 Final remarks 



The approach followed here relies on the relation between the discrete two-particle 
spectrum on the torus and the S-matrix of the corresponding infinite- volume theory 
I10j . Since this relation is only guaranteed to hold up to exponential corrections in 
the volume for massive field theories, we were led to consider massive vector bosons. 
However, in the final QCD matrix element (|2.23p . the dispersion relation of the vector 
bosons only appears to play a secondary role. It would therefore be desirable to have 
a more direct approach to the question of finite-size corrections of the ttq — > 77 decay 
amplitude, presumably using some of the mathematical techniques developed in [9j[2lJ 

In isospin-symmetric QCD, the 77 is perfectly stable, and it would be interesting to 
study the case 77 — > 77 along the same lines we followed here. One must however take 
into account that at -y/i = M v , 77 —> 3ir is kinematically allowed and of the same order 
in the fine structure constant as the 77 pole contribution to light-by-light scattering. 
Unless one makes use of the specific dynamics of pions [23], one is then dealing with 
a case involving an open inelastic channel. The theoretical framework to study such 
problems is under construction [241125] . 
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A Reduction to a one-channel problem 

We consider the elastic scattering of two vector bosons with same helicity. In the final 
state, the vector bosons can have either both positive, or both negative helicity, and 
there are thus two open channels. We follow the notation of Sharpe and Hansen [22], 
The S-matrix is written as 

i(M) B = K((S) B -l), K = ^f- (A.l) 

In the basis B = (\v + v+), \v-V-)), we saw that the scattering amplitude takes the form 

>w»-i- M M -%)■ «-) 

On the other hand, in order to relate A4 to a phase shift, we realize that in the basis 
B' = (\{vv) + ), \(vv)-}), the S-matrix is given by 

(S)B' = diag(e 2i<5< \0). (A.3) 
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Performing the change of basis 
i(M) B =K 



1 1 
-1 1 



D 2iS q 

1 



1 -1 
1 1 



1 
1 



(A.4) 



the scattering amplitude in the basis B thus takes the form (|A.2p with 

e 2iS _ -y 



M = K- 



2i 



(A.5) 



This is precisely the expression given in Eq. (|2.6p . 
B The fermion four-point function 

Using the plane- wave expansion of the quantum fields if) and ijj, we can write the two- 
fermion correlation function as 



^(t,x')i; sl {t',y') il> a (p,-p)$p(p,p) 



(B.l) 



t>t'>o 1 



P p",p' s,s',s" 



F F —Ht' F F 
a P "s" a p's"' e a ps a (-p)s' 



J(p"- X '+p'-y>) 



^JIEpiilEpi 

Now we insert a complete set of states in the remaining matrix element, the contribution 
of a bound state \B) is 



(0\a^ sl ,a^ slll \B)e- Est ' (B\a§ s af_ p y\0). 
The individual matrix elements are given by 



(B.2) 



(0\a n s na p i s ,n\B) 



1 



p" s p. 

,F „F 



V2L* 
1 



:*(p') V'(-p') (<5s"(-)<5s'"(+) + ^s"(+)^'"(-)) ) (B.3) 
(B\ a ; s af_ p)s ,|0) = -jL= (-**(-p)) (*.'(+)*-(-) + <*✓(-)<*.(+)) • (B-4) 
Finally, we obtain the expression for t' — > oo 

(B.5) 



^(t,x')rp s ,(t',y') ^ a (0,-p)^(0,p) 

— Est' 

' / v (~)(p)4 + )(-p ) + v (+)(p)4-)(- ?3 ) 



E e " } z ip " {x '- v,) n-p") ^-\p")u^\-p") + 4+)(p")4" ) (-p // )) • 
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Equivalently, we can transform all fermion fields to the time-momentum representation 
and then obtain 



^ 7 (i,fc)^'(* , ,-<Z / ) M0,~P)M^P)) (B.6) 



e 



--Est' 



e -E k (t-t') _ 



2Ek 



#(-*) (4-)(fc)4^ ) (-fe)+4 +) (fe)4~ ) (- fe )) . 
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